Excited hadrons on the lattice: Baryons 
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We present results for masses of excited baryons from a quenched calculation with Chirally Im- 
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I. INTRODUCTION 

The reproduction of the hadron mass spectrum from 
first principles is an important challenge for lattice QCD. 
Ground state spectroscopy on the lattice is by now a well 
understood problem and impressive agreement with ex- 
periments has been achieved. However, the lattice study 
of excited states ^-^I?| is not as advanced. The reason 
for this is twofold: Firstly, the masses of excited states 
have to be extracted from subleading exponentials in the 
spectral decomposition of two-point functions. Secondly, 
the construction of hadron interpolators which have a 
good overlap with the wave functions of excited states is 
much more challenging than for the ground state. 

Concerning the first issue, the extraction of the signal 
from the subleading exponential, several approaches such 
as constrained fitting or the maximum entropy method 
can be found in the literature [l^, [l3 ■ Here we apply 
the variational method [ij, [lj| j where not only a single 
correlator is analyzed, but a matrix of correlators is used. 
This matrix is built from several different interpolators, 
all with the quantum numbers of the desired state. The 
variational method also incorporates in a natural way 
a solution to the second issue, the wave function of the 
excited states: One uses a set of basis interpolators which 
is large enough to span ground and excited states and 
the variational method finds the optimal combinations of 
them. In principle, no prior knowledge of or assumption 
about the composition of the physical hadron state has 
to be used. 

However, the variational method can succeed only if 
the basis set of hadron interpolators is rich enough to 
span ground and excited states. On the other hand, the 
basis should also be constructed such that it can be im- 
plemented numerically in an efficient way without the 
need for many different quark sources. In this article we 
use a twofold strategy for building our basis interpola- 
tors: We use interpolators with different Dirac structures 
and furthermore compose them using different types of 
smearing for the individual quarks. In particular, wc 
apply different amounts of Jacobi smearing |15J | and in 
this way create "narrow" and "wide" sources. A com- 



bination of these allows for spatial wave functions with 
nodes, which are essential for a good overlap with excited 
states. 

Following a first test of the outlined strategy Pi and 
an analysis of mesons with our method \iU^ . in this pa- 
per we present in detail the results obtained for baryons. 
In the next section we collect the basic equations for the 
implementation of the variational method, detail the con- 
struction of our sources and give an overview of the pa- 
rameters of our numerical simulation. Subsequently we 
discuss effective mass plots, the eigenmodes of the cor- 
relation matrix, as well as the baryon masses and their 
chiral extrapolations. The paper closes with a summary 
and an outlook. 



II. OUTLINE OF THE CALCULATION 
A. The variational method 

As already stated, we use the variational method [ij, 
Il4| to extract the masses of ground and excited states. 
Starting from a set of basis operators Oi, i = l,2...iV, 
we compute the correlation matrix 

Cy(t) = (a.(i)o,(o)). (1) 

In Hilbert space these correlators have the decomposition 



C^At) = J2(0\O,\n){n\O]\0), 



'tM„ 



(2) 



Using the factorization of the amplitudes one can show 
[Tjl that the eigenvalues A^*^) (t) of the generalized eigen- 
value problem 

C(i)w^''^ = X'^''\t)C{to)v^'''> , (3) 

behave as 

A('=)(t) = ^-{t-to)M,^^^(^(^^^(t^to)AM,^^^ (4) 

where Mk is the mass of the k-th. state and AMk is the 
difference to the masses of neighboring states. In Eq. Q 



the eigenvalue problem is normalized with respect to a 
timeslice to < t- 

At this point we remark, that the variational method 
can be generalized to include also ghost contributions as 
they appear in a quenched or partially quenched calcu- 
lation. The fact that ghost contributions play a role also 
for baryons was first stressed in yj. In the spectral de- 
composition Q ghosts appear with a modified time de- 
pendence, possibly including also a negative sign. In [I3 
it was shown that the ghost contribution couples to an 
individual eigenvalue (up to the correction term) in the 
same way as a "proper" physical state. Thus, ghost con- 
tributions are disentangled from the desired states and 
need not be modeled in the further analysis of the expo- 
nential decay of the eigenvalues. 

Let us finally stress that also the eigenvectors of the 
generalized eigenvalue problem ^ contain interesting in- 
formation. If one plots the entries of the eigenvector as a 
function of t, one finds that they are essentially constant 
in the same range of i- values where plateaus of the effec- 
tive mass are seen (compare Fig. EJ- These plateaus can 
be used to optimize the interval for fitting the eigenval- 
ues. Furthermore, the eigenvectors encode the informa- 
tion which linear combinations of the basis interpolators 
couple to which eigenvalue and thus provide one with 
a "fingerprint" of the corresponding states. Comparing 
these fingerprints for different values of the quark mass 
is an important cross-check for the correct identification 
of the states. 
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TABLE I: Dirac structures used for nucleon, E, H and A- 
octet, according to Eqs. (|SJ-@. 



A: 



n-. 



Ap = eabcUa{ulCj^Uc) 



^M = '^abcSaislC"/f,Sc) ■ 



(10) 



(11) 



Here we used vector/matrix notation for the Dirac in- 

(i) (i) 

dices. The different possible choices for Tl and r2 are 
listed in Table n 

Our interpolator for the A (fl) still has overlap with 
both spin-i and spin-|. Thus, we need a projection to 
definite angular momentum. We use the continuum for- 
mulation of a spin- 1 projection for a Rarita-Schwinger 
field: 



1 
3p' 



^^{l ■Plt^Pu+PtJ.1'^1 ■P)i 



B. Dirac structure of the baryon interpolators 

For the baryons we analyze, we use the following op- 
erators with different Dirac structures: 

• Nucleon: 

N^^ = eabcT^\a{ulvfd, ~ dlvfu,) . (5) 



• E: 



SW = ea,j:fua{ulTfs,~slTfu,). (6) 



where p^ is the 4-monicntuni, in our case given by (0, m). 
For each component of the projected A (SI) we compute 
the correlator and average these 2-point functions over 
/i, i^ = 1,2,3. 

Finally, our baryon correlators are projected to definite 
parity using the projection operator P^ = i(l±74). We 
obtain two matrices of correlators: 

-.-tE+ 



zi,^' 



(^T-t)E- 



where we have projected with P+ and 

C-.{t) = -7-.-iE- _7+,-{T-t)E+ 



-Z^,e 



^t,^' 



(12) 



(13) 



when using P . These two matrices are combined to 



• A-octet: 



(7) 



(8) 



A-singlet: 



-(i)„, /jTt.(1) 



T-r(i), 



Al = tabcV'-^'Ua^diV^^'s.-siT^U,) 

■\- cyclic permutations of u, d, s 



(9) 



C{t) 



1 



(C+(t)-C-(T-t)) 



(14) 



to improve statistics. This gives rise to the final corre- 
lator C(t) which we then use in the variational method. 
The positive parity states are obtained from the correla- 
tor at small t running forward in time, while the negative 
parity states arc found at large time arguments, propa- 
gating backward in time with T — t. As expected, the 
correlation matrices C(t) are real and symmetric within 
error bars and we therefore symmetrize the matrices by 
replacing C.,j{t) with C.,j{t) = [Ci.j{t) + Cjt{t)]/2 before 
diagonalization. 



C. Quark sources 

In addition to the different Dirac structures, we con- 
struct the interpolators hsted in the last section from 
quarks with sources created by different amount of smear- 
ing. In particular, we use Jacobi smearing I1!tI | with two 
different sets of parameters (number of smearing steps, 
amplitude of hopping term) to create narrow and wide 
sources. The shapes of these sources are approximately 
gaussian with a ~ 0.27 fm for the narrow source and 
a ^ 0.41 fm for the wide source. Details of the source 
preparation and plots of the source shapes can be found 
in II [13. 

Each quark in our baryon interpolators can cither be 
narrow (n) or wide {w), giving rise to the following eight 
combinations for the sources: 



n{nn) , n(nw) , n{wn) , w{nn) , 
n(ww) , w{nw) , w{wn) , w{ww) 



(15) 



In this notation the order of the quark fields is under- 
stood as in Eqs. (O - Hll|) and the parentheses indicate 
which quarks are combined in the diquark combination. 
Since the smearing used here is a purely scalar operation, 
the assignment of quantum numbers, as given in the last 
subsection, remains unchanged. 

Taking into account the different Dirac structures dis- 
cussed in the last section, we can work with 24 different 
interpolators for nucleon, E, S and A-octet. For the A- 
singlet and the A we have only one Dirac structure and 
consequently a total of 8 different interpolators. We re- 
mark that in the final analysis not all interpolators are 
used. We prune the maximal correlation matrix and re- 
move some of the correlators that couple only weakly to 
the physical states or add no new information, thus en- 
hancing the numerical noise. The criterion for the selec- 
tion of the interpolators is the optimization of the quality 
of the plateaus in the effective mass 
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(16) 



D. Parameters of the simulation 

We work with quenched gauge coiifigurations gener- 
ated with the Liischer-Weisz action ^^. We use two 
sets of lattices, 20^^ x 32 and IG'^ x 32, at couplings 
P = 8.15 and /3 = 7.90 corresponding to lattice spac- 
ings of a = 0.119 fm and a = 0.148 fm, determined from 
the Sommer parameter in [T^ . Thus for both lattices we 
have a spatial extent of i ~ 2.4 fm. The two different val- 
ues of the lattice constant a allow us to assess the cutoff 
dependence. The parameters of the gauge configurations 
are collected in Table HTl 

Our quark propagators are computed using the Chi- 
rally Improved (CI) Dirac operator [l^. The CI oper- 
ator is a systematic approximation of a solution of the 
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confs. 


a[fm] 


a-^lMeV] 


20'' X 32 
16^x32 


8.15 
7.90 


100 
100 


0.119 
0.148 


1680 
1350 



TABLE II: Parameters of our simulation. We list the lattice 
size, the inverse coupling /3 the number of configurations, the 
lattice spacing a and the cutoff a~^ . 



Ginsparg- Wilson equation J23 with good chiral prop- 
erties 21]. We work with several quark masses in the 
range am — 0.02, ... 0.2, leading to pion masses down to 
350 MeV. For setting the strange quark mass we use the 
K-meson with the light quark mass extrapolated to the 
chiral limit. 

Our quark sources are placed at i = and the general- 
ized eigenvalue problem Q is normalized at to — lo- The 
final results for the baryon masses were obtained from a 
fully correlated fit to the eigenvalues. The errors we show 
are statistical errors determined with single elimination 
jackknife. 



III. RESULTS 
A. Effective masses, eigenvectors and fit ranges 

Let us begin our presentation with a discussion of ef- 
fective masses (|16ll for the nucleon system. For posi- 
tive parity the combination of the 6 operators n{wwp^' , 
w{wnp^' , w{ww)^^' , n{wwY^' , w{wny^' , w{w'w)^^' , (the 
upper index denotes the choice of Dirac structures ac- 
cording to Table P) gives the strongest signal. For neg- 
ative parity we used the 4x4 correlation matrix built 
from n{nnp^\ w{nnp^\ n{nny'^\ w{nnp'^' . 

In Fig. n we compare the effective mass plots for pos- 
itive and negative parity nucleons from our two lat- 
tices at different values of the bare quark mass; am — 
0.05,0.1,0.2 for 16^ x 32 and am = 0.04,0.08,0.16 for 
20'^ x 32. These numbers were chosen such that they 
give rise to approximately equal pion masses for the two 
lattice spacings used. The plots also contain the nucleon 
masses in lattice units as obtained from a correlated fit of 
the propagator (horizontal bars giving the central value 
plus and minus the statistical error). The figure shows 
clear long plateaus for the ground state masses, while 
the signals for excited states have larger error bars and 
shorter plateaus. Furthermore the quality of the data 
decreases as the quarks become lighter - a feature well 
known in lattice spectroscopy. 

Another important piece of information comes from 
the eigenvectors. In Fig. [3 we show the 6 entries of the 
lowest three eigenvectors corresponding to ground, first 
and second excited state (top to bottom) in the positive 
parity nucleon channel. Again we compare the results for 
our two lattice sizes using quark mass values which give 
rise to essentially the same pion mass. For each value of t 
the respective eigenvectors are normalized to unit length. 
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FIG. 1: Effective mass plots for nucleon ground and excited states. We compare the results from our coarse (16"^ x 32, a = 0.148 
fm, aruq = 0.05,0.1,0.2 top to bottom), and fine (20^ x 32, a = 0.119 fm, arug = 0.04,0.08,0.16) lattices. The solid lines are 
the results from correlated fits of the eigenvalues. They represent the fit result plus and minus the corresponding error. 



It is interesting to note that the eigenvectors are only 
weakly dependent on t (actually this can be shown from 
the generalized eigenvalue problem). The entries form 
plateaus which are very long for the ground states but 
also for the excited states often contain 4 to 8 values of 
t. Typically these plateaus extend at least over the same 
number of i-values as the effective mass plateaus - often 
they are even longer by one or two points. 

As in the case of effective masses, the formation of the 
eigenvector plateaus indicates that the channel is domi- 
nated by a single state. Thus, the eigenvector plateaus 
provide an important tool for the reliable identification of 
the i-intervals where the eigenvalues can be used for a fit. 
Indeed, sometimes it is the eigenvectors which prevent 
one from fitting "quasiplateaus" in the effective mass. 
Due to relatively large statistical errors in the effective 
masses, the data sometimes resemble a plateau and it 
is only the absence of a plateau in the corresponding 
eigenvectors which allows us to conclude that a quasi- 
plateau is not conclusive. We implemented this strategy 
and now fit the eigenvalues only where we see also eigen- 
vector plateaus. 

We finally remark that the values for the eigenvectors 
are almost exactly the same for the two values of the 



cutoff we consider (the left hand side plots are for a — 
0.148 fm, the right hand side is for a = 0.119). Although 
the entries of the eigenvectors cannot be expected to scale 
(they are linear combinations of matrix elements of our 
interpolators with the physical states), it is reassuring for 
the application of the method that no large discrepancies 
are observed. 



B. Nucleon 

As already discussed in the last section, the pos- 
itive parity nucleon masses were extracted from the 
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6x6 correlation matrix of n{ww) 
n{'wwp^', w{wnp^' , 'w{wwp^> , while for negative parity 
the 4x4 correlation matrix built from n{nnY^' , w{nnY^' , 
n{nnp'^' , w{nnp'^' was used. Of course these combina- 
tions were used for all quark masses. For the positive 
parity ground state we could determine the mass for all 
our quark masses. For the excited nucleon states of pos- 
itive parity the combined assessment of effective masses 
and eigenvector plateaus did not allow for a trustworthy 
extraction of the corresponding nucleon masses for the 
two smallest quark masses. 
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FIG. 2: Eigenvectors for nucleon ground and excited positive 
parity states. From top to bottom we show the eigenvector 
components of ground, first and second excited state. 



We identify two excited states of positive parity which 
have not very different masses for the whole quark mass 
region where we see a signal. This is consistent with our 
previous observation on a smaller lattice ^]. These arc 
two physically distinct states since they are observed in 
different eigenvalues of the correlation matrix and the 
corresponding eigenvectors are orthogonal. Some addi- 
tional efforts are required to properly identify the nature 
of our quenched excited states. We follow the strategy 
of Ref. [3, i.e., we trace the states from the heavy quark 
region towards the physical limit. 

In the heavy quark region, where we obtain the best 
signals, the quenching and chiral symmetry effects are 
less important and the naive quark picture is adequate. 
Then we know a-priori, that there must be two approxi- 
mately degenerate excited states of positive parity. The 
first one is a member of the 56-plet of SU(6) and the 
second one belongs to the 70-plet. In the excited 56-plet 
state, as well as in the ground state 56-plet (the nucleon) , 
all possible quark pairs have positive parity. Then it fol- 
lows that the signal from the nucleon, as well as from the 
excited 56-plet state, can be seen with those interpola- 
tors that contain two-quark subsystems of positive parity 
(these are the ones with i — 1,3 from the Table I). On 
the other hand, the positive parity 70-plet state contains 
both positive and negative parity two-quark subsystems, 
and can be seen with the i = 2 interpolator, where the 
"diquark" has negative parity. This picture is confirmed 
in the heavy quark limit of our results. If we construct 



our correlation matrix with the i — 1 and/or i = 3 in- 
terpolators, we find both the ground state (the nucleon) 
and two excited states of positive parity, while only one 
state is observed with the i = 2 interpolator. This state 
corresponds to the positive parity excited state. 

Using the "fingerprint" from the eigenvectors, we are 
able to trace these states from the heavy quark region, 
where their physical nature can be safely identified, to the 
light quark region (down to m-^ = 450 MeV) , where they 
still remain approximately degenerate. Clearly these sig- 
nals, extrapolated to the physical region, remain essen- 
tially higher than the experimental states N(1440) and 
N(1710) (cf. Fig.|2Jl. Phenomenologically the latter states 
are ascribed to the 56-plet and 70-plet, respectively. 

The discrepancy between our results and the exper- 
imental numbers is probably partly due to quenching, 
where a significant part of chiral physics is absent. Also 
finite volume effects cannot be excluded (our physical 
volume is 2.4 fm and large finite volume effects can be 
anticipated for excited states (g)- 

Note that the perturbative gluon exchange between va- 
lence quarks, characteristic of the naive constituent quark 
model, is adequately represented in the quenched calcula- 
tion. The discrepancy of our results with the experimen- 
tal ones hints that it is the chiral physics, partly missing 
in quenched QCD, that could shift both positive p arity 
excited states (and especially the Roper state) down J22 . 

Our results for the nucleons are presented in Fig. |31 
The left plot is for positive parity, the right for negative 
parity. Filled symbols are used for the 16'^ x 32, a = 0.148 
fm lattice, open symbols for 20^ x 32, a = 0.119 fm. The 
filled circles represent the experimental masses. 

The results for the positive parity ground state (left 
plot, downward pointing triangles) agree well with the 
experimental value (for the chiral extrapolation of our 
data see Subsection F). Furthermore, the data show al- 
most no cutoff effects. For the first excited state (circles) 
the results for the two values of the cutoff differ by about 
one sigma, while for the second excited state (upward 
pointing triangles) the two data sets agree. However, 
both excited states extrapolate to values about 20-30% 
larger than the experimental numbers. 

For negative parity we mainly fit ground and first ex- 
cited states. Only for the two largest quark masses on 
the finer lattice we can extract the second excited mass. 
We find that the lowest two states are nearly degenerate, 
but extrapolate to the physical masses within error bars 
(compare Subsection F). Cutoff effects are clearly seen 
only for small quark masses. Since the negative parity 
ground and first excited state are nearly degenerate, we 
checked that they are indeed different by inspecting the 
eigenvectors and following their behavior down from the 
heavy quark region. Entries of the eigenvectors at quark 
mass am — 0.06 are shown for our 20^ x 32 lattice in 
Fig. 21 In contrast to the positive parity excited states, 
the negative parity states fit the experimental data well. 
This is expected since the negative parity states have the 
mixed flavor-spin symmetry [21]i?5 and experience only 



positive parity 



negative parity 



O 2-0 






N(1710) 
; N(1440) 



^' 



• N(938) 
] , 



3:5: 



^ 



0.80 

(mj' [GeV]' 



5 
5 



l" 


a = 0. 


l" 


a = 0. 


i'-' 


a = 0. 


i" 


a = 0. 


i" 


a = 0. 


i'-' 


a = 0. 



14Sfin ■ 
14Sfin 
148fin 
119fin ■ 
119fin 
119fin 
I 



1.20 



N(2090) 



• N(1650) 

• N(I535) 



« 



l" 


a = 0. 


x" 


a = 0. 


x" 


a = 0. 


x" 


a = 0. 


x" 


a = 0. 



0.80 

(mj' [GeV]' 



14Sfm 
14Sfin 
119fm . 
119fin 
119fin ' 



1.20 



.4 

.0 

.6 
.4 
2 
.0 
.8 



FIG. 3; Ground and excited state nucleon masses versus M^ for our two lattices. Filled symbols are used for 16'^ x 32, a — 0.148 
fm, open symbols for 20 x 32, a = 0.119 fm. The left hand side plot shows the positive parity states, the right hand side is for 
negative parity. The experimental data are included as filled circles. 



small chiral effects Ba. 



C. E and H 

Those E and S resonances which belong to the octet 
are structurally identical to the nucleon: only one and 
two, respectively, of the light quarks are replaced by a 
strange quark. Consequently, their analysis and also the 
results are only a variation of what has been found for 
the nucleon system. We use the same combination of 
interpolators in the 6x6 (for positive parity) and 4x4 
(negative parity) correlation matrices as we did for the 
nucleons. 

We present our results for the octet E and S masses in 
Fig. El As for the nucleon system, the positive parity E 
and S ground states are compatible with the experimen- 
tal numbers and essentially no cutoff effects are visible. 
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FIG. 4: Eigenvectors for nucleon ground and first excited 
negative parity states. The data are for our 20'^ x 32 lattice 
at am = 0.06. 



Concerning the excited positive parity states, only the 
first excited states show notable cutoff effects, while the 
masses of the second excited states from the two lattices 
are compatible within error bars. For the S, where at 
least the first excited state is classified, our data extrap- 
olate to a number which is about 20 % larger than the 
experimental result, similar to the nucleon case. 

For negative parity, we find two nearly degenerate 
states which show clear cutoff effects for the smaller 
quark masses. For the E the data are compatible with 
the known states. For the negative parity S our data 
extrapolate to two states near 1800 MeV (see also the 
discussion in Subsection F). 



D. A 

For A we have considered two different kinds of inter- 
polators; one which is a pure flavor singlet and one which 
has mainly overlap with a flavor octet. 

For the flavor octet A we obtain results which are sim- 
ilar to the results of the other flavor octet baryons. Even 
the same combination of sources used for N, E and S 
turns out to be the optimal one also for the A octet chan- 
nel. 

For the flavor singlet A we are mainly interested in the 
ground state in both parity channels. We have therefore 
used only a single interpolator, the one where all quarks 
are smeared narrowly (choosing a different smearing com- 
bination does not change the results). 

The interesting observation is that while the negative 
parity flavor-singlet state extrapolates to the mass which 
is essentially higher than A(1405), which is consistent 
with previous quenched lattice results, the flavor-octet 
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FIG. 5: Same as Fig. now for E and 



negative parity ground state signal is consistent with the 
A(1405). Within the simple quark model picture the 
negative parity pair 1/2", A(1405) - 3/2", A(1520) is a 
flavor-singlet. However, starting from the early Dalitz' 
work it is also understood that at least a significant part 
of A(1405) could be due to KN physics (2^]. The KN 
bound state system can couple to the flavor-octet interpo- 
lator and our results hint at the KN nature of A(1405). It 
would be very interesting to study the 3/2~, A(1520) res- 
onance and to see whether it is a flavor-singlet or flavor- 
octet state. 



E. A2 and Q^ 

As already discussed in the previous section, our in- 
terpolators for A 2 have to be spin projected to obtain 



correlators of states with definite quantum numbers. Af- 
ter the projection we are left with a set of 8 interpola- 
tors which differ only in the smearing combination of the 
quarks. From these we have chosen different subsets and 
found that the dependence on the chosen subset is only 
marginal. In the end, we decided to use the combina- 
tions n{nn), w{nn), n{n'w), w{nw), n{ww), wlww) for 
both parity channels. 

In Fig. [7| we present the results for the A 2 and il^ 
masses. The positive parity states are shown in the left 
plot, the right plot is for negative parity. The vertical 
lines in both plots mark the values of m^ corresponding 
to the physical strange quark mass, which has been de- 
termined from a fit to the K-meson mass in a separate 
calculation on the two lattices. At these values of the 
pion mass we extract the masses for the il^ resonance 
from our results for the A2. It is remarkable that the 
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FIG. 6: Ground and excited state masses obtained from our A octet (upper plots) and A singlet (lower plots) interpolators. 



ground state f^s lies right on top of the experimental 
value. 

The results for the positive parity ground states of A 2 
show significant discrepancies with the experimental re- 
sults. However, this is not unexpected and has already 
been observed by other groups Q- The Ropcr-like state, 
A(1600), is not reproduced either. In both cases the most 
probable explanation would be a lack of the proper chiral 
dynamics in quenched QCD. Given the fact that the il 
ground state is perfectly reproduced, one concludes that 
this missing chiral dynamics becomes especially impor- 
tant at the quark masses below the strange quark mass. 

On the negative parity side we could only reliably fit 
the ground state and only on the fine lattice do our data 
reach the strange quark mass such that the mass of the 
negative parity fl^ can be determined. Extrapolation to 
the physical limit is consistent with A(1700). 



F. Chiral extrapolations for the fine lattice 

Where the data are sufficient, we perform a chiral ex- 
trapolation of our results. For excited states the form of 
the chiral extrapolation is not known from chiral pertur- 
bation theory and we extrapolate linearly in m^. Since 
in this paper the focus is on the excited states, the ex- 
trapolation for the ground states is also kept simple - we 
use second order polynomials in jn^ there, which is the 
structure of the leading terms in quenched chiral pertur- 
bation theory [2j. Since for some of the states we still 
observe cutoff effects, we extrapolated only the data from 
the finer lattice. 

For positive parity the results of the chiral extrapola- 
tion are presented in the left plot of Fig. |S1 We remark, 
that the numbers for the fi are obtained by an interpola- 
tion to the strange quark mass. While the ground states 
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FIG. 7; A masses versus A/^. The vertical lines mark the values of M^ corresponding to the physical strange quark mass. 



come out reasonably well for a quenched calculation, the 
results for the excited states are systematically 20 % - 
25% above the experimental numbers (where known). 
The most likely explanation is that quenching removes 
some important piece of chiral physics, which is actually 
responsible for the proper mass of excited positive parity 
states. Significant finite volume effects cannot be ruled 
out either. 

For negative parity states (right plot of Fig. ^ the 
results are compatible with the experimental numbers 
(where known), although the statistical errors are larger. 
Also here we cannot exclude that cutoff effects push our 
numbers up a little bit, but from the comparison of the 
results on our two lattices we estimate that this effect 
is not larger than the statistical error. Again the result 
for Q is obtained from an interpolation to the strange 
quark mass. One may expect that quenching effects are 
essentially smaller for the negative parity channel states 
than for the positive parity excited states. This is ex- 
pected a-priori, since all low-lying negative parity states 
have the mixed flavor-spin symmetry [21]fs and hence 
are affected by the chiral dynamics only slightly (except 
for the A(1405) ) 0. 



IV. SUMMARY 

In this article we presented a quenched spectroscopy 
calculation of excited baryons using the variational 
method. We use interpolators with different Dirac struc- 
tures. Furthermore each quark can either have a narrow 
or a wide source such that the states can have nodes in 
their spatial wave function. 

For the positive parity baryons we find that the ground 
state masses are compatible with the experimental num- 
bers, while for the excited states the masses are system- 



state 



Q, positive parity, first excited state 
Q, negative parity, ground state 
E, negative parity, ground state 
H, negative parity, first excited state 



Mass [MeV] 



2300(70) 
1970(90) 
1780(90) 
1780(110) 



TABLE III: Collection of our final results for some states not 
classified by the Particle Data Group [23] . 



atically 20% - 25% above the experimental numbers. We 
believe that the failure to reproduce the masses of the 
positive parity excited baryons is indeed mainly due to 
quenching where a significant part of chiral physics is 
missing. Large finite volume effects cannot be excluded 
either. 

For negative parity, we find that our masses are in rea- 
sonable agreement with the experimental numbers, al- 
though here our statistical errors are larger and a further 
lowering of our results for a lattice with a higher cutoff 
cannot be excluded. 

In some of the channels we analyze, the corresponding 
baryons are not yet classified [23 ■ For four of these chan- 
nels we believe that our data are strong enough to quote 
the final results as a prediction: The first excited positive 
parity O state, the negative parity Q ground state, and 
the ground and first excited negative parity S states. 

The two fl states are included in this list since at the 
strange quark mass the chiral dynamics is less important 
and also our results do not need to be extrapolated to the 
chiral limit. Concerning the two negative parity S states 
we believe that the good results of the structurally very 
similar negative parity nucleons and S baryons justify the 
prediction of the mass of the negative parity ground and 
first excited state in the S channel. Our final numbers 
for the masses of the four states are listed in Table IIIII 



10 









positive parity 






2.5 


- 










- 


2.0 


- 5 


I 


5 


5 


i 


I 


J 1.5 


- ™ 


IMii 






o 


--: 




- 


^^ 


O 






- 


1.0 


- n 


z 


r^ 


O 

A 


A^ 


Q." ; 



negative parity 




FIG. 8: Chiral extrapolation of our results. The left plot is for positive parity, the right for negative parity. The horizontal bars 
represent the experimental numbers (where known), indicating also the error. For our results we use circles for ground states, 
squares and diamonds for the first and second excited states respectively. The shaded square symbol for the excited A represents 
the chiral limit of the data from the singlet interpolator. Filled symbols are used for those states where no corresponding state 
is listed in the particle data book. 



Ackno'wledginents 

C.G. acknowledges interesting discussions on the vari- 
ational method with Rainer Sommer and Peter Weisz. 
The calculations were performed on the Hitachi SR8000 
at the Leibniz Rechenzentrum in Munich and the JUMP 
Cluster at the Zentralinstitut fiir angewandte Mathe- 
matik in Jiihch. We thank the LRZ and ZAM staff 



for training and support. L.Y.G. is supported by 
"Fonds zur Forderung der Wissenschaftlichcn Forschung 
in Osterreich", FWF, project P16823-N08. This work is 
supported by DFG and BMBF. C.G. and C.B.L. thank 
the Institute for Nuclear Theory at the University of 
Washington for its hospitality and the Department of 
Energy for partial support during the completion of this 
work. 



[2] 
[3] 



[5] 



[1] N. Mathur et ai, Phys. Lett. B 605, 137 (2005); F. X. 

Lee et ai, Nucl. Phys. Proc. Suppl. 119, 296 (2003). [8] 

S. Sasaki, Prog. Theor. Phys. Suppl. 151, 143 (2003). 
S. Sasaki, T. Blum and S. Ohta, Phys. Rev. D65, 074503 
(2002). 

[4] D. Brommel et al. [Bern-Graz-Regensburg CoUabora- [9] 

tion], Phys. Rev. D 69, 094513 (2004); Nucl. Phys. Proc. 
Suppl. 129, 251 (2004). [10] 

T. Burch et al. [Bern-Graz-Regensburg Collaboration], 
Phys. Rev. D 70, 054502 (2004); Nucl. Phys. Proc. Suppl. [11] 
140, 284 (2005); Nucl. Phys. A 755, 481 (2005); PoS [12] 
LAT2005, 075 (2005). 

[6] K. Sasaki, S. Sasaki and T. Hatsuda, Phys. Lett. B 623, [13] 

208 (2005); K. Sasaki and S. Sasaki, PoS LAT2005, 060 [14] 
(2005); S. Sasaki, K. Sasaki, T. Hatsuda and M. Asakawa, [15] 
Nucl. Phys. Proc. Suppl. 129, 212 (2004); Nucl. Phys. 
Proc. Suppl. 119, 302 (2003); [16] 

J. M. Zanotti et al, Phys. Rev. D 68, 054506 (2003); 
W. Melnitchouk et al, Phys. Rev. D 67, 114506 (2003); [17] 

Nucl. Phys. Proc. Suppl. 109 A, 96 (2002); F. X. Lee et 
al, Nucl. Phys. Proc. Suppl. 106, 248 (2002); F. X. Lee 
and D. B. Leinweber, Nucl. Phys. Proc. Suppl. 73, 258 



[7] 



(1999). 

K. J. Juge et al , "arXiv:hep-lat/0601029' S. Basak et al 

[Lattice Hadron Physics Collaboration (LHPC)], Phys. 

Rev. D 72, 074501 (2005); Phys. Rev. D 72, 094506 

(2005); arXiv:hep-lat/0601034 

D. Guadagnoli, M. Papinutto and S. Simula, Phys. Lett. 

B 604, 74 (2004). 

T. Burch et al [Bern-Graz-Regensburg Collaboration], 

arXiv:hep-lat/0601026 PoS LAT2005, 097 (2005). 

G.P. Lepage, Nucl. Phys. (Proc. Suppl.) 106, 12 (2002). 

M. Asakawa, T. Hatsuda and Y. Nakahara, Prog. Part. 

Nucl. Phys. 46, 459 (2001). 

C. Michael, Nucl. Phys. B 259, 58 (1985). 

M. Liischer and U. Wolff, Nucl. Phys. B 339, 222 (1990). 

S. Giisken, et al, Phys. Lett. B 227, 266 (1989); C. Best 

et al, Phys. Rev. D 56, 2743 (1997). 

T. Burch, C. Gattringer, L.Y. Glozman, C. Hagen and 

C.B. Lang, Phys. Rev. D 73, 017502 (2006). 

M. Liischer and P. Weisz, Commun. Math. Phys. 97, 59 

(1985), Err.: 98, 433 (1985); G. Curci, P. Menotti and 

G. Paffuti, Phys. Lett. B 130, 205 (1983), Err.: 135, 516 

(1984). 



11 



[18] C. Gattringer, R. Hoffmann and S. Schaefer, Phys. Rev. 
D 65, 094503 (2002). 

[19] C. Gattringer, Phys. Rev. D 63, 114501 (2001); C. Gat- 
tringer, I. Hip and C. B. Lang, Nucl. Phys. B 597, 451 
(2001). 

[20] P.H. Ginsparg and K.G. Wilson, Phys. Rev. D 25, 2649 
(1982). 

[21] C. Gattringer et al. [Bern-Graz-Regensburg Collabora- 
tion], Nucl. Phys. B 677, 3 (2004); Nucl. Phys. Proc. 
Suppl. 119, 796 (2003); C. Gattringer, Nucl. Phys. Proc. 
Suppl. 119, 122 (2003). 

[22] L. Ya. Glozman and D. O. Riska, Phys. Rep. 268, 263 



(1996); L. Ya. Glozman, W. Plessas, K. Varga, R. Wagen- 

brunn, Phys. Rev. D 58, 094030 (1998); L. Ya. Glozman, 

Nucl. Phys. A 663, 103 (2000). 
[23] For a review, see R. H. Dalitz (p. 748), in Review of 

Particle Physics, (D. E. Groom et al), Eur. Phys. J. C 

15, 1 (2000). 
[24] J.N. Labrenz and S.R. Sharpe, Phys. Rev. D 54, 4595 

(1996). 
[25] S. Eidclman et al. (Review of Particle Physics), Phys. 

Lett. B 592, 1 (2004). 



